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Abstract 

We present a theoretical approach which allows to determine the number 
and orbital character of the conducting channels in metallic atomic contacts. 
We show how the conducting channels arise from the atomic orbitals having a 
significant contribution to the bands around the Fermi level. Our theory pre- 
dicts that the number of conducting channels with non negligible transmission 
is 3 for Al and 5 for Nb one-atom contacts, in agreement with recent exper- 
iments. These results are shown to be robust with respect to disorder. The 
experimental values of the channels transmissions lie within the calculated 
distributions. 
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Metallic contacts of atomic dimensions, that can be produced by means of scanning tun- 
neling microscope and break-junctions techniques, have been the object of great attention 
in the last few years [IJ. According to the scattering approach, the electronic transport in 
these mesoscopic structures can be described in terms of independent "conducting chan- 
nels" characterized by certain transmission coefficients which vary between zero and one 
||. The complete understanding of their transport properties requires the knowledge of 
the transmission coefficients along each conducting channel. This information, which is not 
accessible by usual conductance measurements, has for the first time been obtained in recent 
experiments on superconducting Al contacts ||. The possibility of extracting the individ- 
ual transmissions from measuring the superconducting I — V characteristic is based on the 
extreme sensitivity of the subgap structure to small changes in these parameters 

A remarkable conclusion of this experimental study is that, although the total conduc- 
tance of an Al contact on the first plateau (pressumably a one-atom contact ||) can be 
close to one quanta of conductance, this situation does not correspond to a nearly open 
single channel, but rather to a situation with three partially open channels. The question 
naturally arises on the microscopic origin of this phenomena for the specific case of Al. More 
generally, one would like to know what is the number of conducting channels for a given 
metal in a given contact geometry. The aim of this letter is to provide theoretical insight 
into these questions. 

Having a system of atomic dimensions it seems natural to choose an atomic orbital basis 
for the description of its electronic structure. This choice has proven useful in the context of 
STM theory 0. Furthermore, the use of a local basis in combination with Green function 
techniques provides an efficient way for obtaining the transport properties in terms 

of microscopic parameters. In an atomic orbital basis the electronic Hamiltonian adopts the 
usual tight-binding form 

H — ^ia c i aa C iaa + J]] ti a jfjC ia o .Cjfj tCr , (1) 
ia,ir ia^jf3,cr 

where i, j run over the atomic sites and a, f3 denote the different atomic orbitals (the number 
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of orbitals at each site will be denoted by N or b). The hopping elements Uajp are assumed 
to connect first-neighbor sites only. There exist in the literature several parametrization 
procedures for determining the tight-binding Hamitonian [|9|,|T0|| which are known to accu- 



rately reproduce the band structure of bulk materials. As a starting point we follow the 
parametrization proposed in Ref. |TD| considering as a minimal basis for each metal those 



atomic orbitals having a significant contribution to the electronic DOS around the Fermi 
energy, Ep. Thus, for the case of superconducting metals of groups III and IV, like Al, Pb, 
etc, only s and p orbitals need to be considered; while for transition metals like Nb d orbitals 
have to be included. 

In an atomic contact the local environment in the neck region is very different to that 
of the bulk material and therefore the use of bulk parameters in the Hamiltonian requires 
some justification. In the first place, the inhomogeneity of the contact geometry can produce 
large deviations from the approximate local charge neutrality that typical metallic elements 
must exhibit. Within the tight-binding approximation this effect can be corrected imposing 
local charge neutrality through a self-consistent variation of the diagonal parameters €i a 



|IIfl. As discussed below, this self-consistency in the neck region turns out to be crucial 
for the correct determination of the conducting channels. Regarding the hopping elements, 
tiajfi, although we shall initially consider them as being equal to the bulk values in order to 
represent a neck geometry with bulk interatomic distances; we shall show that the results 
are robust with respect to fluctuations in the hopping elements induced by disorder in the 
atomic positions. 

An idealized geometry for a one- atom contact is depicted in Fig.l. It consists of a close- 
packed fee structure grown along the (111) direction (hereafter denoted as z direction), 



starting from a central atom [|12|]. This structure is connected to two semi- infinite crystals 
describing the metallic leads, N being the number of atomic layers within the neck region. 
By taking different values of N we can describe both the cases of long and short necks. 

In order to obtain the d.c. current for a constant bias voltage applied between the 
leads, it is most convenient to use non-equilibrium Green function techniques [|T3 1 . Within 



this formalism the current can be written with an expression formally equivalent to that of 
scattering theory |TJ| by treating the coupling between the central region and the leads as a 
perturbation jnj . The current between the left lead and the central region is then given by 



/oo 
T(E,V)[f L (E)-f R (E)]dE, (2) 
-oo 



2e r°° 
~h 



where fL )R are the Fermi-distribution functions for the (left, right) leads and T(E, V) is an 
energy and voltage dependent transmission probability, which can be written in terms of the 
matrix elements of the (retarded, advanced) Green function matrix G r ' a (E) = [E±iO + —H]~ 1 



as 



T(E, V) = 4Tr 



lm± L (E-e^) x 



G r m (E)lmE R (E + e^)G a N1 (E) 

In this expression G r 1N and G a Nl are matrices whose elements are the Green functions 
connecting layers 1 and N, 1£l,r being self-energy matrices describing the coupling of the 
central region to the leads. These matrices have a dimension equal to the number of bonds 
connecting the central region to the leads (Mx^) and have a simple expression in terms of 
the Green functions of the uncoupled leads, gtajp- 



ia.jB 

J k,l£L,R;y,6 

The giajpS can be evaluated numerically by standard decimation techniques |T^| . 

The voltage range which is probed in the experiments with superconducting contacts 
is of the order of the gap parameter A (typically a few tenths of meV) 0. For this bias 
range normal metallic systems will behave ohmically. Even when the atomic-size contacts 



exhibit resonances around Ep [Ojl, their width will in general be much larger than A and 
linearization of Eq. (2) is appropriate. In this linear regime the contact normal conductance 
can be expressed as G = (2e 2 /h) T(Ep, 0). By using the cyclic property of the trace, G can 
be written in the Landauer form G = (2e 2 /h) Tr[t(Ep)t^ (Ep)], where 



1/2 



t(E) = 2 lmZ L (E) G r 1N (E) ImZ R (E) . (3) 



1/2 



The existence of (ImS) 1 / 2 as a real matrix is warranted by ImS being positive definite. 
Moreover, W is an hermitian matrix having Ml real eigenvalues, Tj, which are bounded by 
zero and one |LjJ. Associated with these eigenvalues there will be Ml eigenvectors, which 
in our model are linear combinations of the atomic orbitals in the layer which is in contact 
with the left lead. These eigenvectors define the way in which the atomic orbitals contribute 
to each conducting channel. 

Although the dimension of ift can be arbitrarily large depending on the size of the central 
region, the actual number of conducting channels (those with a non-vanishing transmission 
eigenvalue Tj) are limited by the number of orbitals in the narrowest section of the neck 
(N or b when having a single atom contact). This fact can be shown by the following simple 
argument. As the division between "central region" and leads is somewhat arbitrary, one 
could always redefine the leads for the geometry of Fig. 1 in such a way that the new central 
region would only consist of the central atom. Then, the new lead self-energy matrices R 
would have a dimension of just N or b and the new transmission matrix would only admit 
N or b eigenmodes. Current conservation along each conducting channel ensures that the 
nonvanishing eigenvalues Tj and T[ must be the same. 

The above simple argument already allows an estimate of the maximum number of 
relevant conducting channels in a one- atom contact. Thus, for an sp-like metal like Al, this 
number should be typically four, while for a transition metal like Nb (having a negligible 
weight of p orbitals at Ep) this number would be of order six. As discussed below, this rough 
estimate should be taken as an upper bound. The actual number of conducting channels can 
be smaller as some of the channels can carry no current due to symmetry considerations. 

Let us first consider the case of an Al one- atom contact, which is the one analyzed 
experimentally in Ref. 0. Al contacts have also been addressed theoretically in [18|]. The 
atomic configuration for Al, 3s 2 3p 1 , gives rise to a conducting sp band with three electrons 
per atom. The bulk DOS at Ep has important contributions from both 3s and 3p orbitals, 



the 3s level being located ~ 7eV below the 3p level [10]. While the above simple argument 
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would predict a maximum number of four channels for a one-atom contact, the self-consistent 
calculation for the ideal geometry yields only three channels with nonvanishing transmission. 
In Fig. 2 the transmission eigenvalues of the ideal (111) contact are shown as a function of 
energy for the two extreme cases of a short (N = 1) or a long (N — > oo) neck geometry. 
There are certain features that are common to both cases: i) there are three channels having 
a significant transmission around Ep, the fourth one being closed for every energy; ii) the 
total transmission is close to one around Ep increasing to almost three at higher energies; 
iii) there is a non-degenerate mode which is widely open for almost every energy; iv) the 
second transmission eigenvalue is two- fold degenerate and has a small value around Ep. 

The channels can be classified according to the orbital character of the eigenvectors on the 
central atom. In this way, the non-degenerate channel, widely open at Ep, is a combination 
with an amplitude ~ 0.63 on the s and ~ 0.77 on the p orbitals along the z direction. 
The one with zero transmission corresponds to the orthogonal combination of these two 
orbitals. On the other hand, the two degenerate modes correspond to combinations of p 
orbitals on a plane perpendicular to the z direction. While the symmetry properties of the 
neck geometry are responsible for the decoupling between the s — p z and p x — p y orbitals, 
the approximate fulfillment of the condition lmT,' ss lmE' pzPz = (lmS' spz ^ accounts for the 
presence of a non-conducting channel (details will be given elsewhere). 

The orbital character and energy dependence of the transmission eigenvalues is similar 
for the case of Pb, which is also an sp-like metal. However in this case the Fermi level lies 
in the region where the three channels are almost open, giving rise to a total transmission 
larger than 2. 

So far the possible effect of disorder in the atomic positions has been disregarded. We 
have studied this effect by introducing random fluctuations in the atomic positions of the 
idealized structure, assuming the distance dependence on the hopping parameters as sug- 
gested by Harrison |9|]. Although there are certain features like the two-fold degeneracy 
which, as expected, disappear with the inclusion of disorder, the gross features found for 
the ideal geometry are nevertheless robust. This fact is illustrated in the histograms shown 
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in Fig. 3, which demonstrate that the decomposition of the total transmission consists of a 
widely open channel with Tj between 0.6 and 0.9 and two low transmissive channels with 7] 
between 0.1 and 0.3. These predictions are consistent with the experimental results for the 
first conductance plateau [[|. The fourth channel always has an extremely small transmission 

m < io- 4 ). 

We have also analyzed the conducting channels of a Nb one-atom contact, as an example 
of transition metal contacts. As commented above, in this case the maximum number of 
conducting channels is expected to be six due to the fact that the DOS around Ep mainly 
arises from the contributions of 5s and 4d orbitals [|Kj]. The idealized one-atom contact 
geometry yields in this case a total transmission between 2 and 3 depending on the number 
of layers in the central region. The channel decomposition shows that this total transmission 
is mainly built up from the contribution of five conducting channels. As can be observed 
in Fig. 4, the d bands cause a strong energy dependence of the transmission eigenvalues, 
with typical energy scales of the order of 0.5eV. The s and d z 2 orbitals hybridize strongly 
and give rise to the conducting channel with the highest transmission around Ep (mode 1 
in Fig. 4). The almost closed channel (mode 6 in Fig. 4) corresponds to the orthogonal 
combination of these two orbitals. There also appears a two-fold degenerate channel with 
transmission ~ 0.7 and another two-fold degenerate channel with transmission ~ 0.3. Both 
the values of the total transmission and the number of relevant conducting channels are in 



good agreement with preliminary experimental results [|i~9f| . 

Our theory and the experimental results show that atomic contacts of metals which have 
an important contribution from p and d orbitals, do not necessarily, even in the one-atom 
case, exhibit an integer number of perfectly transmitting modes. This situation is at variance 
with that of simple metals like Na, Au, Ag, etc, which can be described by a single s like 
band around Ep and exhibit well defined quantized conductance steps, at least at the lower 
plateaus ||20|| . Within our theory, in the one-atom contact geometry, simple metals would 



have a single conducting channel. The transmission of that channel is strongly pinned at 
one due to the charge neutrality condition |TT[| . 



In conclusion, we have presented a theoretical analysis of the conducting channels in 
metallic atomic-size contacts. We have shown that the number and character of these chan- 
nels are determined by the orbital electronic structure and the local atomic environment 
around the neck region. For the case of sp-like metals like Al and Pb we predict the pres- 
ence of 3 conducting channels in a one- atom contact, in good agreement with the available 
experimental data |§. For one-atom contacts of a transition metal like Nb we expect the 
presence of 5 conducting channels due to the contribution of d orbitals. This result has been 



confirmed by recent experiments [IS 
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FIGURES 

FIG. 1. Idealized geometry for a one-atom contact. The layers are numbered from 1 to N 
starting from the left lead. 

FIG. 2. Transmission eigenvalues as a function of energy for Al one-atom contacts in the two 
extreme cases of short (a) and long (b) necks. The solid curve corresponds to the non-degenerate 
sp z mode and the dotted curve corresponds to the two-fold degenerate p x — p y mode. 

FIG. 3. Typical distributions of the transmission eigenvalues for an Al one-atom contact when 
disorder in the atomic positions is included (the maximum fluctuation in the hopping parameters 
is of the order of 100% with respect to the bulk values, which corresponds to variations of the order 
of 20% in the interatomic distances). The two modes corresponding to the two-fold degenerate 
eigenvalue in the idealized case exhibit a similar distribution. 

FIG. 4. Transmission eigenvalues as a function of energy for a Nb one-atom contact for the 
long neck case. Note the two-fold degeneracy of modes 2-3 and of 4-5. 
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